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Abstract 

We discuss novel QCD phenomena in pA collisions at LHC energies such as 
the possibility to invesigate the behaviour of hard processes at x as low as 
10~^, signals for onset of the black body limit in the hard processes and in the 
nucleon fragmentation, multijet production. 

1. Physics motivations 

The proton-ion colhsions at LHC will be qualitatively different from those at fixed target energies [] In 
the soft interactions two prime effects are a strong increase of the mean number of "wounded" nucleons 
- from about 6 to about 20 in pA central collisions (due to increase of ainei{NN)) and a 50 % increase 
of the average impact parameters in pp collisions (as manifested in the shrinkage with increase with 
energy of the diffractive peak in the elastic pp collisions) making it possible for a proton to interact 
often simultaneously with nucleons which have essentially the same longitudinal coordinate but different 
impact parameters. At the same time geometry of pA collisions at LHC should be more close to that in 
the classical mechanics because increase of cross sections combined with the suppression of the inelastic 
diffraction in pp scattering, leads to the suppression of fluctuations in the value of effective NN cross 
section (and hence the suppression of inelastic shadowing corrections in pA scattering), to a very strong 
suppression of inelastic diffraction in pA collisions as compared to the fixed target FNAL energies [^. 
In the hard interactions the prime effect is a strong increase of the gluon densities at small x both in 
projectile and ion. Gluon densities can be measured at LHC up to ultra small x a 10^^^'' where PQCD 
will be probed in a new domain [|l]]. As a consequence of blackening of the interaction of proton partons 
with nuclear partons at small xa one should expect disappearance in the proton fragmentation region 
of hadrons with small transverse momenta, and hence a strong suppression of nonperturbative QCD 
effects Thus the high transverse momentum proton fragmentation region is a natural place to search 
for unusual QCD phases which should live long time because of the Lorentz dilatation of time. In spite 
of a small coupling constant, gluon densities of heavy ions become larger in a wide range of x and 
and impact parameters than that permitted by the conservation of probability within the leading twist 
approximation [ |T2| , |l|, ^, ^, ^. Thus a decrease of the strong coupling constant with virtuality would be 
insufficient to support applicability of PQCD to the hard small x processes in the kinematical domain 
where gluon fields may appear strong enough for nonperturbative QCD vacuum to become unstable and 
where equations of QCD may need to be modified. The possibility for novel QCD regime and therefore 
new phenomena are maximized in this QCD regime. f\ At the same time interpretation of these novel 
effects would be much more definitive in pA collisions than in the heavy ion collisions. 

New classes of strong-interaction phenomena which occur within the short-distance, near-the- 
light-cone space-time region and which could be associated with perturbative QCD or with interface of 
perturbative and nonperturbative QCD are much more probable in pA collisions than in pp collisions. 
The densities of partons, and of energy-momentum, may be high enough to "burn away" the nonpertur- 
bative QCD vacuum structure in a cylinder of a radius ~ Ifm and length ~ 2Ra, leaving behind only 

'An extensive discussion of the novel pA physics at LHC was presented in |jl|]. Here we summarize few key points of this 
analysis and extend it to consider effects of the small x high field regime. 

^A similar phenomenon of the suppression of small transverse momenta for another quantity, the nuclear light cone wave 
function at small x, is discussed in the saturation models 1^ ^• 

'it is important to distinguish between parton densities which are defined within the leading twist approximation only |^] 
and may increase with energy forever and physical quantities - cross sections/structure functions - whose increase with energy 
at given impact parameter is restricted by the probability conservation. 



a dense partonic fluid governed by new QCD dynamics with small but effective short-distance large 
couplings. If there is a large enhancement of hard-collision, gluon-induced processes, then there should 
be an enhancement of hard multiparton collisions, heavy-flavor production [|l]]. The study of the x, pi 
and A-dependence of multiparton interactions, of the forward charm and beauty hadron production (in 
the direction of the fragmenting proton) should be especially incisive to address these questions. 

Estimates of the kinematical region for the new QCD dynamics in proton (nucleus)-heavy nucleus 
collisions based on the LO formulae for the dipole cross section and the conservation of probability 
indicate that for the central impact parameters for A ~ 200 and 2bGeV^ for quarks xa < 

5 • 10^^ and for gluons xa < 10~^. 



2, Measurements of nuclear parton densities at ultra small x 

One of the fundamental issues in high energy QCD is the dynamics of the hard interactions in the small x 
kinematics. Depending on the resolution scale (Q^) one investigates here either the leading twist effects 
or the regime of strong color fields. Since the cross-sections for hard processes increase roughly linearly 
with A (except for very small x and relatively small Q"^ where the counting rates are high anyway), even 
short runs with nuclear beams will produce data samples sufficient to measure parton distributions with 
statistical accuracy better than 1 % practically down to the smallest x which is kinematically allowed: 
X ~ 10~^ for quarks and x ~ 10~^ for gluons [^. The systematic errors for the ratios of nuclear and 
nucleon structure functions are also expected to be small since most of these errors cancel in the ratios 
of the cross-sections, cf 

There are several processes which could be used to probe the small x dynamics in pA collisions. 
They include the Drell-Yan pair production, study of dijets, jet -i- photon, charm production, diffractive 
exclusive production of three jets [p3|], multijet events. It is worth emphasizing that for these measure- 
ments it is necessary to detect jets, leptons and photons at rapidities close enough to the nucleon rapidity. 
In this y-range, the accompanying soft hadron multiplicities are relatively small, leading to a soft parti- 
cle background comparable to that in pp scattering. In fact, for \ymax — y\ < (2 — 4) the background 
level is likely to be significantly smaller than in pp collisions due to suppression of the leading particle 
production in pA scattering (for a thorough discussion of these reactions see [|l|]). 

The measurements of the dimuon production appear to be more feasible experimentally for the 
forward angles than studies of the hadron production due to small energy losses for muons. In the 
kinematics where the leading twist dominates one can study both the quark small x distributions (from 
the cross section at relatively small transverse momenta: p1{/j,'^fi~) <^ M^(//+^~), and the gluon 
distribution (from the measurements of the cross section at pt{^~^ ^ M (fi~^fi~) [[T6|]). Note also that 
for the forward region a dimuon background from the charm production is likely to be rather small since 
muons in average carry only 1/3 of the energy in the charm decays. 



3. Signals for onset of the black body limit 

A fast increase with energy of structure functions of heavy nuclei contradicts to the probability conser- 
vation for a wide range of parton momenta and virtualities for the LHC energy range. On the contrary 
the rapid, power-like growth of structure function of a nucleon observed at HERA may continue forever 
because of the diffuse edge of a nucleon [||]. At present new QCD dynamics of the regime of high parton 
densities where multiparton interactions are inhibited remains a subject of lively debates. PreQCD black 



body regime for the structure functions of heavy nuclei [IC], QCD black body regime for hard processes 
off heavy nuclei [^], saturation of parton densities at given impact parameter [^, color glass condensate 
regime [Q], turbulence and related scaling laws are possible options. For the estimates of the pattern 
of the novel phenomena we shall use in this text the formulae of black body regime because they are 
generic enough for the expectations of the onset of new QCD dynamics within the existing theoretical 



approaches. 

The use of nuclei because of their uniform density and a possibiUty to select scattering at the 
central impact parameters allows one to enhance the high parton density effects. Indeed, in the small x 
regime a hard collision of a parton of the projectile nucleon with a parton of the nucleus occurs coherently 
with all the nucleons at a given impact parameter. The coherence length l/2mNX by far exceeds 
the nuclear size: in the kinematic regime accessible at LHC it can reach up to 10^ fm (in the nucleus 
rest frame). In the rest frame of the nucleus this can be visualized in terms of the propagation of a 
parton in high density gluon fields over much larger distances than is possible with free nucleons. In 
the Breit frame it corresponds to the fact that small x partons cannot be localized longitudinally to 
anything close to the Lorentz-contracted thickness of the nucleus. Thus low x partons from different 
nucleons overlap spatially, creating much larger parton densities than in the free nucleon case. This 
leads to a large amplification of the multiparton hard collisions expected in QCD at small x; see e.g. 
[||, H]. Eikonal approximation [11], constraints from the probability conservation [12, IJ] indicate that 
new QCD dynamics should reveal itself at significantly larger x than in the proton and in more striking 
effects. Recently extensive studies of hard diffraction were performed at HERA, for review see [|l^]. This 
allows to evaluate taming due to the leading twist nuclear shadowing which turns out to be large both 
in the gluon and quark channels and significant up to very large Q^, - for the recent discussion see [^. 
However the leading twist taming being large is still insufficient to restore the probability conservation 
leaving a lot of room for the violation of leading twist approximation. It maybe more serious than merely 
due to enhancement of the high twist effects as the whole picture of interactions can qualitatively change. 
Really, the common wisdom is that smallness of guarantees that bare particle is free within a hadron 
but this remarkable property which confirmed experimentally in many hard processes is violated in small 
X phenomena. One may think about this phenomenon in term of an effective violation of the asymptotic 
freedom in this limit. It is expected to occur in heavy nuclei at x at least ten times larger than in the 
proton. Hence in the LHC kinematics the region where novel QCD phenomena are important will extend 
at least two orders of magnitude in xa for a large range of virtualities. 

The space time picture of the black body limit in the proton-nucleus collisions can be seen best in 
the rest frame for the nucleus. A parton belonging to the proton emits a hard gluon (virtual photon) long 
before the target and interacts with the target in a black regime releasing the fluctuation, e.g. a Drell-Yan 
pair. This leads to a qualitative change in the picture of the proton-nucleus interactions for the partons 
with Xp,pt satisfying the condition that 

XA = (1) 

XpSNN 

is in the black body kinematics for the resolution scale pt < p'1'^''"{xa)- Here pj'^ ' (x^) is maximum pt 
for which the black body approximation is applicable. In the kinematics of LHC « 4(pj '' ' )^ can be 
estimated by using formulae derived in [12]. At minimal xa it may appear as large as (p^'^' ^ C^^a))^ = 
15GeV^ . All the partons with such Xp will obtain pt{jet) ~ pj'^ ' (xyi) leading to the multijet pro- 
duction. The black body regime will extend down in Xp with increase of the incident energy. For LHC 
for Pt < 3GeV/c this regime may cover the whole region of Xp > 0.01 where of the order ten par- 
tons resides. Hence in this limit most of the final states will correspond to multiparton collisions. For 
Pt ^ 2GeV/c the region extends to Xp > 0.001. At LHC collisions of such partons correspond to the 
central rapidities. Dynamics of conversion of the high pt partons with similar rapidities to hadrons is 
certainly a collective effect which requires a special consideration. 

Inclusive observables For the calculation of the total cross sections the logic similar to the one 
used in the consideration of 7* - nucleus scattering [^ should be applicable. In particular for the case of 

''in the black regime highly virtual probe interacts at the same time with several high kt partons from the wave function of 
a hadron |^. This is qualitatively different from the asymptotic freedom limit where interaction with only one highly virtual 
parton are important. Hence the interacting parton is far from being free when probed by a hard probe of a large coherence 
length = l/2mNX. 



the total cross section of the dimuon production can be evaluated as: 

-F2p(xp, • -^M'' ■ 2T^R\hi{xo/xA), (2) 



dxAdxp 9 M2 / g^2- 

for large but not too large dimuon masses, M^. In particular we estimate M'^{bbl){xA = 10~^ 
GOGeV^. Here A'-factor has the same origin as in the leading twist case, but it should be smaller since 
it originates from the gluon emissions only from the parton belonging to the proton, Ra is the nuclear 
radius, and xq is maximal xq for which the black body limit is valid. Hence the expected Af ^ depen- 
dence of the dimuon production is qualitatively different in this case than in the case of pp scattering 
where scattering at large impact parameters may mask the black body limit contribution. Also the xa 
dependence in this limit becomes pretty weak. 

The study of the pt distribution of the dimuons may provide also a signal for the onset of the black 
body regime. Similar to the case of pt distribution of leading partons in the deep inelastic scattering 
one should expect a broadening of the pt distribution of the dimuons in the black body limit as compared 
to the DGLAP expectations, see [ [l7| ] for a calculation of this effect in the color glass condensate model. 

With further decrease of xa formulae of black limit will probably overestimate cross section be- 
cause the interaction with heavy nucleus of sea quarks and gluons in the proton would become black 
also. 

The onset of the black body limit will lead also to gross changes in the hadron production: much 
stronger drop with xp of the spectrum in the proton fragmentation region accompanied by a gross pt 
broadening of the spectrum as well as the enhancement of the hadron production at smaller rapidities, 
see discussion in section I 



It is worth emphasizing here that onset of the black body limit for ultra small xa will pose limits 
on using the forward kinematics for the measurements of the parton densities at larger xa- Indeed as soon 
as one tries to use large enough Xp for the nuclear parton density measurements one would have to take 
into account that collisions of these partons with nuclear partons are always accompanied by a significant 
Pt broadening associated with the interaction with "black component" of the nucleus wave function. As 
a minimum this would lead to a significant pt broadening of the pt distribution of the produced hard 
system. One can also question the validity of the leading twist expansion for the cross sections integrated 
over the transverse momentum of the produced system. To investigate these phenomena an ability to do 
measurements at fixed xa and different Xp would be very important. 



4. Mapping of the three dimensional nucleon parton structure 

The systematic studies of hard inclusive processes during the last two decades have led to a pretty good 
understanding of the single parton densities in nucleons. However very little is known about multi- 
parton correlations in nucleons which can provide critical new insights into the dynamics of the strong 
interactions, and allow to discriminate between different models of nucleons. Such correlations may be 
generated, for example, by the fluctuations of the transverse size of the color field in the nucleon leading, 
via color screening, to correlated fluctuations of the densities of gluons and quarks. 

A related source of correlations is QCD evolution, since a selection of a parton with a given 
X, may lead to a local (in transverse plane) enhancement of the parton density at different x values. 
Also, practically nothing is known about possible correlations between the transverse size of a particular 
configuration in the nucleon and the longitudinal distribution of partons in this configuration. 



4.1 Multi-jet production and double parton distributions 



It was recognized already more than two decades ago [18 1 that the increase of parton densities at small x 
leads to a strong increase of the probability of nucleon-nucleon collisions where two or more partons of 
each projectile experience pair-vice independent hard interactions. Although the production of multijets 



through the double parton scattering mechanism was investigated in several experiments Ill9|, gO|] at pp, pp 
colliders, the interpretation of the data was hampered by the need to model both the longitudinal and the 
transverse partonic correlations at the same time. The studies of proton-nucleus collisions at LHC will 
provide a feasible opportunity to study separately the longitudinal and transverse partonic correlations in 
the nucleon as well as to check the validity of the underlying picture of multiple collisions. 

The simplest case of a multiparton process is the double parton collision. Since the momentum 
scale pt of a hard interaction corresponds to much smaller transverse distances 1/pt in the coordinate 
space than the hadronic radius, in a double parton collision the two interaction regions are well separated 
in the transverse space. Also in the cm. frame pairs of partons from the colliding hadrons are located in 
pancakes of thickness < {l/xi + l/x2)/pc.m.- So two hard collisions occur practically simultaneously 
as soon as xi,X2 are not too small and hence there is no cross talk between two hard collisions. A 
consequence is that the different parton processes add incoherently in the cross section. The double 
parton scattering cross section, being proportional to the square of the elementary parton-parton cross 
section, is therefore characterized by a scale factor with dimension of the inverse of a length squared. The 
dimensional quantity is provided by the nonperturbative input to the process, namely by the multiparton 
distributions. In fact, because of the localization of the interactions in transverse space, the two pairs of 
colliding partons are aligned, in such a way that the transverse distance between the interacting partons of 
the target hadron is practically the same as the transverse distance between the partons of the projectile. 
The double parton distribution is therefore a function of two momentum fractions and of their transverse 
distance, and it can be written as T{x,x',b). Actually F depends also on the virtualities of the partons, 
Q^, Q'^, though to make the expressions more compact we will not write explicitly this dependence. 
Hence the double parton scattering cross section for the two "two —>■ two" parton processes a and /3 in 
an inelastic interaction between hadrons a and b can be written as: 

m I ^ , ,,, .',2, 



a£){a,P) = — J Ta{xi,X2;b)aa{xi,x[)ai3{x2,x'2)Tb{x'i,x'2;b)dxidx'idx2dx'2d b (3) 

where m = 1 for indistinguishable parton processes and m = 2 for distinguishable parton processes. 
Note that though the factorization approximation of Eq.(||) is generally accepted in the analyses of the 
multijet processes and appears natural based on the geometry of the process no formal proof exists in the 
literature. As we will show below the study of the A-dependence of this process will allow to perform a 
stringent test of this approximation. 

Here to simplify the discussion we neglect small non-additive effects in the parton densities, which 
is a reasonable approximation for 0.02 < x < 0.5. In this case we have to take into account only b- space 
correlations of partons in individual nucleons. 

One has therefore two different contributions to the double parton scattering cross section: = 
aj-, + ajj. The first one, erf, interaction with two partons of the same nucleon in the nucleus, is the 
same as for the nucleon target (the only difference being the enhancement of the parton flux) and the 
corresponding cross section is 

cj]j = aD j d^BT{B) = Aao, (4) 



where 



/oo r 
dzpA{r),jT{B)d^B = A (5) 



is the nuclear thickness, as a function of the impact parameter of the hadron-nucleus collision B. 

The contribution to the term in rA{x'i,x'2, b) due to the partons originated from different nucleons 
of the target, af^, can be calculated solely from the geometry of the process by observing that the nuclear 
density does not change within a transverse scale (6) ^ Ra- It rapidly increases with A oc / T^{B)d'^B. 
Using information from the CDF double scattering experiment[^] on the mean transverse separation of 



partons in a nucleon, one finds tliat tlie contribution of the second term should dominate in the case of 



proton - heavy nucleus collisions: crjj/aj-, « .68(j4/12)-^^ [21|. Hence one expects a stronger than 
oc A increase of the multijet production in pA collisions at LHC. Measurements with a set of nuclei 
would to measure the double parton distributions in nucleons and also to check the validity of the QCD 
factorization for such processes which appears natural but which so far was not derived in pQCD. An 
important application of the discussed process would be to investigate transverse correlations between 
the nuclear partons in the shadowing region. This would require a selection of both partons from the 
nucleus in the shadowing region, xa < ~ 10~^. 

As we already mentioned in section |] partons with sufficiently large Xp satisfying Eq.|T] are ex- 
pected to interact with a probability of the order one with small x partons of nucleus. As a result a hard 
collision of the partons with sufficiently large Xp and xa > 0.01 will be accompanied by production of 
a minijet/minijets at the black body kinematics with pt ~ p^'^ ' (a;/i). Most of sufficiently fast partons of 
the protons will generate minijets leading to a strong suppression of cross section of the events with only 
two jets. 

Detectors with sufficiently large rapidity acceptance would allow to detect events originating from 
the triple parton collisions at xa > Xsh and large enough pt where pt broadening related to the blackbody 
limit effects are sufficiently small. This provide stringent tests of the dynamics of the hard interactions 
and provide additional information on two parton correlations as well as unique information on triple 
parton correlations. 

Other opportunities with multi-jets include 

• Probing correlations between partons in the nucleus at high densities, i.e. for xi a,X2 a 10~^, 
which would provide qualitatively new information about the dynamics of nuclear shadowing and 
the presence of possible new condensates of partons. 



• Studying the accompanying soft hadron production (cf. the discussion in Section |4.2| ), which 
would allow measurement of the transverse size of a proton configuration containing partons 
xi,X2- In particular, the production of 4 jets with two jets at large x in the proton fragmentation 
region xi + X2 > 0.5 may provide another way to look for point-like configurations in nucleons, 
see [ p3| ] and [|l|] for discussion of some other options. 

4.2 Proton-ion collisions probe transverse nucleon structure 



As we discussed in []15[], one supposes that in some subclass of events the distribution of constituents 
in the initial proton may be unusually local in the transverse (impact parameter) plane when the proton 
collides with the ion. If this is so, its effective cross-section per nucleon will be greatly reduced, perhaps 
all the way to the perturbative-QCD level. If the effective cross-section of such a point-like configuration 
goes below 20mb, there will be an appreciable probability that it can penetrate through the center of a 
Pb ion and survive. This would lead to a highly enhanced yield of diffractive production of the products 
of the point-like configuration in collisions with heavy ions. 

Not only might the properties of the final-state collision products depend upon the nature of the 
transverse structure of the proton primary at arrival at the collision point, but even the conventional parton 
distributions may also be affected. For example, let us consider, a nucleon as a quark and small diquark 
connected by a narrow QCD flux-tube. It should be clear that if the flux-tube is at right angles to the 
collision axis at arrival, then the valence partons will have comparable longitudinal momentum or x. On 
the other hand, if the flux tube is parallel to the direction of motion, then one of the valence systems will 
have very large x, and the other very small. This happens because in this case the internal longitudinal 
momenta of the valence systems, in the rest frame of the projectile proton, are in opposite directions. 
Therefore the smallness of the configuration is correlated with the joint j;-distribution of its constituents. 

'The A-dependence of the ratio of o-ijax in the kinematics where only one of the nuclear partons has xa < x^h is 
practically the same as for the case when both nuclear partons have x > Xsh- 



This kind of nonfactorization may be determined by the study of the perturbative-QCD processes of 
dilepton, direct-photon, or dijet production as a function of the centrality of the coUision, multiplicity of 
soft hadrons, as well as a function of the atomic number. |^ Naturally such studies would require also 
experimental investigation of the production of the soft hadrons as a function of impact parameter at 
LHC energies where it may differ quite strongly from the one observed at fixed target energies. 

One possible kinematics where a strong correlation is expected is when a parton with large x{x > 
0.6) is selected in the proton. The presence of such a parton requires three quarks to exchange rather large 
momenta. Hence one may expect that these configurations have a smaller transverse size and therefore 
interact with the target with a smaller effective cross-section (7e//(x). Suggestions for such a dependence 
of the size on x are widely discussed in the literature. Using as a guide a geometric (eikonal type) picture 
of pA interactions and neglecting (for simplicity) shadowing effects for nuclear parton densities one can 
estimate the number of wounded nucleons u{x,A) in events with a hard trigger (Drell-Yan pair, 7-jet, 
dijet,...) as a function of (jg// [|^]: 

u{x, A) = l + a,ff{x)^^ J T\h) d\ (6) 



where the nuclear density per unit area T{h) is defined in Eq.( 11). At LHC for average inelastic pPb 
collisions and a^jf ~ CTineiiPP) Eq.§ leads to 1/ 10. This is somewhat larger than the average number 
of wounded nucleons in pA collisions, due to the selection of more central impact parameters in events 
with a hard trigger. 

A decrease of the effective cross-section for large x, say, by a factor of 2, would result in a com- 
parable drop of the number of particles produced at central rapidities as well as in a smaller number of 
nucleons produced in the fragmentation of the nucleus, in the nucleus fragmentation region. 

4.3 ^-dependence of the particle production 

Particle production in the proton fragmentation region The j4-dependence of hadron production in 
the proton fragmentation region remains one of the least understood aspects of hadron-nucleus interac- 
tions. Practically all available data are inclusive and correspond to the energies where the cross section 
of inelastic NN scattering is about three times smaller than at LHC. They indicate that the cross-section 
is dominated by the production of leading particles at large impact parameters, where the projectile in- 
teracts with only one or two nucleons of the nucleus. As a result very little information is available about 
hadron production at the central impact parameters which are most crucial for the study, for example of 
AA collisions. Theoretical predictions for this region are also rather uncertain. 

In eikonal type models, where the energy is split between several soft interactions, one may expect 
a very strong decrease of the yield of the leading particles. The dependence is expected to be exponential 
with path length in the nucleus; their mean energy is attenuated exponentially. On the other hand, if the 
valence partons of the projectile do not lose a significant amount of their initial momentum, as suggested 



by the models motivated by perturbative QCD, see review in []23[], the spectrum of leading particles may 
approach a finite limit for large A and central impact parameters [24]. Indeed in this case the leading 
partons will acquire significant transverse momenta and will not be able to coalescence back into leading 
baryons and mesons as it seemingly happens in the nucleon-nucleon collisions. As a result they will 
fragment practically independently leading to much softer distributions in the longitudinal momentum 



for mesons and especially for baryons [g4j]. 

We discussed above that the strength of the interaction of fast partons in the nucleon (xp > 10^^ 
for LHC) with heavy nuclei at the central impact parameters may approach the black body limit and 

*The impact-parameter dependence is accessible in thie proton - nucleus collisions via a study of nuclear fragmentation into 
a number of channels, for an extensive discussion see [|l|]. 



that under this scenario all these partons will acquire large transverse momenta ~ 



h.b.l. 



xa). Q This 

would lead to a very significant pt broadening of the spectrum of leading hadrons [26]. Probably the 
most feasible way to study this effect will be to measure the pt distribution of the leading neutrons for 
the collisions at the central impact parameters, which will closely follow the pt distribution of leading 
quarks [21 \ It is worth emphasizing that a strong pt broadening which should increase with increase of 
X F clearly distinguishes this mechanism of the suppression of the leading hadron spectrum from the soft 
physics effect of the increase of the total pp cross-section by nearly a factor of three from fixed target 
energies. 



Particle production in the central region The Gribov-Glauber model including inelastic screening 
effects, but neglecting the final state interactions between the hadrons produced in different Pomeron 
exchanges, as well as neglecting the interactions of these hadrons with the nuclear target, leads to a 
prediction that at high energies and rapidities y, such that \yp — y| » 1, the inclusive spectrum of 
produced hadrons should be proportional to that in the inclusive cross-section of hN scattering: 

daP+^-^+^{y,PT) ^ ^^ daP+^^'^+^iy,pT) 
dy d^pT dy d'^pT 

Since ainei{pA) / ainei{pN) oc A^^'^, Eq.( ^ implies that the multiplicity of particles produced in the 
inelastic pA collisions should increase with A, roughly as A^^^. Data at fixed target energies do not 
contradict this relation, but the energy is too low for an unambiguous interpretation. Effects of splitting 
the energy between multiple soft A^A^ interactions are much larger in the nucleus-nucleus collisions than 
in proton-nucleus collisions so checking whether similar formula is valid for AA collisions would require 
energies much higher than those available at RHIC. 

The opposite extreme is to assume that interactions of partons in the black body regime become 
important in a large range of Xp > 10^2 jn this case multiple scatterings between the partons, their 
independent fragmentation together with associated QCD radiation (gluon bremsstrahlung) may produce 
the major part of the total entropy and transverse energy and will lead to enhancement of the production 
of hadrons at central rapidities as compared to the the expectations of the Gribov-Glauber model. 



Particle production in the ion fragmentation region For the rapidities close to the nucleus rapidity, 
there are indications of the contribution of slow hadron re-interactions, which lead to an increase of the 
multiplicity of nucleons with momenta 1.0>pAr>0.3 GeV, as compared with expectations neglecting 



final state re-interactions (a factor of ~ 2 for heavy nuclei) [|28|]. As we have discussed above, it is doubt- 
ful that the Gribov-Glauber picture is correct at LHC energies, so one may expect significant deviations 
from the expectations based on Eq.( 0). It is likely that color quark system with rapidities close to yA will 
be excited along the cylinder of the radius ~ 1 fm and length Ra- This may result in unusual properties 
of the hadron production in the nucleus fragmentation region including production of exotic multiquark, 
multigluon. . . states. 

It must be emphasized that a quantitative understanding of the physics of particle production in 
pA collisions in this phase-space region is a prerequisite for the understanding the corresponding phe- 
nomenology in A A collisions. 

In conclusion, studies of proton-nucleus scattering at LHC will allow to study strong interactions 
in the high field domain over extended region of rapidities. New phenomena will be especially prominent 
in the proton fragmentation region, though they will also extend to the central region and to the nucleus 
fragmentation region. 

'Hence in the black body limit pt broadening of the partons should be much larger than at low energies where pt broadening 
is consistent with the QCD multiple rescattering model |E3||, see [E5|| for the discussion of the matching of these two regimes. 



It is our pleasure to thank J.Bjorken for the joint study of the pA physics for the FELIX project 
and numerous discussions of the forward physics. We also thank GIF and DOE for support. 



References 

[1] E. Lippmaa et al. [FELIX Collaboration], "FELIX: A full acceptance detector at the LHC," CERN- 
LHCC-97-45, J. Phys. G 28 (2002) Rl 17. 

[2] L. Frankfu rt, G. A. Miller and M. Strikman, Phys. Rev. Lett. 71 (1993) 2859 [a rXiv:hep- 
ph/9309285]. 

[3] A. H. Mueller, |arXiv:hep-ph/0111244 . 



[4] E. lancu, A. Leonidov and L. McLerran, |arXiv:hep-ph/020227C . 
[5] R.P. Feynman, Photon - hadron interactions, Benjamin, 1972. 

[6] L. Frankfurt, V. Guzey and M. Strikman, J. Phys. G 27 (2001) R23 [ kXiv:hep-ph/0010248D ; 



L. Frankfurt, V. Guzey, M. McDermott and M. Strikman, arXi v : hep-ph/0 1 04252 



[7] L. Alvero, J. C. CoUins, M. Strikman and J. J. Whitmore, Phys. Rev. D 57 (1998) 4063 [a rXiv:hep- 
ph/97 10490]. 

[8] M. Arneodo, A. Bialas, M.W. Krasny, T. Sloan, M. Strikman In ^Hamburg 1995/96, Future physics 
at HERA* 887-926. 

[9] L. Frankfurt, V. Guzey, M. McDermott and M. Strikman, Phys. Rev. Lett. 87 (2001) 192301 
| |arXiv:hep-ph/0104154| ]. 

[10] V. N.Gribov, Zh. Eksp. Teor. Fiz. 57 (1969) 1306. 

[11] A. H. Mueller and Jian-wei Qiu, Nucl. Phys. B 268 (1986) 427; A. H. Mueller, Nucl. Phys. B 335 
(1990) 115. 

[12] L. Frankfurt, W. Koepf, M. Strikman, Phys. Rev. Phys. Rev. D 54 (1996) 3194; L. L. Frankfurt, 
M. F McDermott and M. Strikman, JHEP 9902 (1999) 002 [ ]hep-ph/98123T6| l. 

[13] H.Abramowicz, L.Frankfurt, M. Strikman Pubhshed in SLAC Summer Inst. 1994:539-574. 

[14] H. Abramowicz and A. Caldwell, Rev. Mod. Phys. 71 (1999) 1275 [ lhep-ex/9903037| ]. 

[15] L. Frankfurt and M. Strikman, "A new hard QCD phenomenon - proton diffraction into three jets". 
Contribution to the workshop "Hard Probes in Heavy Ion Collisions at the LHC", and references 
therein. 



[16] S. Catani et al, "QCD,' arXiv:hep-ph/0005025 



[17] F Gelis and J. Jalilian-Marian, |arXiv:hep-ph/0208141 



[18] C. Goebel, F Halzen and D.M. Scott, Phys. Rev. Dll (1980) 2789; N. Paver and D. Treleani, Nuovo 
Cimento A70, (1982) 215; B. Humpert, Phys. Lett. B131 (1983) 461 (1983); B. Humpert and R. 
Odorico, Phys. Lett. 154B (1985) 211; T. Sjostrand and M. Van Zijl, Phys. Rev. D36 (1987) 2019. 

[19] T. Akesson et al. (AFS Collaboration), Z. Phys. C34 (1987) 163; J. Alitti et al. (UA2 Collaboration), 
Phys. Lett B268 (1991) 145; F Abe et al. (CDF Collaboration), Phys. Rev. D47 (1993) 4857. 

[20] F Abe et al. (CDF Collaboration), Phys.Rev.Lett. 79 584 (1997); Phys. Rev. D56 (1997) 3811. 



[21] M. Strikman and D. Treleani, Phys. Rev. Lett. 88 (2002) 031801 QarXivihep-ph/Oll 1468| ]. 
[22] L. Frankfurt and M. Strikman, Nucl.Phys.250 (1985) 143. 

[23] R. Baier, D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50 (2000) 37 [a rXiv:hep- 
ph/0002198]. 

[24] A. Berera, M. Strikman, W. S. Toothacker, W. D. Walker, J. J. Whitmore, Phys. Lett. B403 
(1997) 1. 

[25] Yu. Dokshitzer, L. Frankfurt, and M. Strikman, in preparation. 

[26] A. Dumitru and J. JaUUan-Marian, Phys. Rev. Lett. 89 (2002) 022301 [ ]arXiv:hep-ph/0204028| | . 
[27] A. Dumitru, L. Gerland, J. Jalilian-Marian, and M. Strikman, in preparation. 
[28] L. Frankfurt and M. Strikman, Phys. Rept. 76 (1981) 215. 



